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^i.^-1 ^ We present an accurate nucleon-nucleon (NN) potential based upon chiral effective 

Lagrangians. Tiie model includes one- and two-pion exchange contributions up to 
chiral order three and contact terms (which represent the short range force) up to 
order four. Within this framework, the NN phase shifts below 300 MeV lab. energy 
and the properties of the deuteron are reproduced with high-precision. This chiral 

C ^ ' NN potential may serve as a reliable starting point for testing the chiral effective 

^Nj ^ field theory approach in exact few-nucleon and microscopic nuclear many-body 

calculations. 

T-H ■ 

> 

f~>. ' 1 Introduction 

vn ' 

^^ ' One of the most fundamental problems of nuclear physics is to derive the force 

[ , between two nucleons from first principles. A great obstacle for the solution 

of this problem has been the fact that the fundamental theory of strong in- 
^—v , teraction, QCD, is nonperturbative in the low-energy regime characteristic 

for nuclear physics. The way out of this dilemma is the effective field theory 
concept which recognizes different energy scales in nature. Below the chiral 
symmetry breaking scale, A^ « 1 GeV, the appropriate degrees of freedom are 
CJ> , pions and nucleons interacting via a force that is governed by the symmetries 

of QCD, particularly, (broken) chiral symmetry. 

The derivation of the nucleariJorcc from chira]|X;ffective field theory was 
pioneered by Weinberg)il, Ord6fipz,B and van KxikikBcl Important contrihiitions 
were made by Robilotta et al.^ Kaiser et aZ.j2rLli3 and Epelbaum et al.n As a 
}J] ' result, efficient methods for deriving the nuclear force from chiral Lagrangians 

5y \ have emerged. Also, the quantitative nature of the chiral NN potential has 

improved.El Nevertheless, even the currently 'best' chiral NN potentials are 
too inaccurate to serve as a reliable input for exact few-nucleon calculations 
or miscroscopic nuclear many-body theory. 
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The time has come to put the chiral approach to a real test in microscopic 
nuclear structure physics. Conclusive results can, however, be produced only 
with a 100% quantitative NN potential based upon chiral Lagrangians. For 
this reason, we have embarked on a program to develop a NN potential that 
is based upon chiral effective field theory and reproduces the NN data with 
about that-apnc. equality as the high-precision NN potentials constructed in 

the iggo'sMm 

In Sees. 2 to 5, we will develop, step by step, the chiral NN potential; and 
in Sec. 6 we will shown that we have achieved our goal. 

2 Effective Chiral Lagrangian 

The effective chiml ttN Lagrangian is given by a series of terms of increasing 
chiral dimension ,E3 

C-„N = C^N + ^^N + ^itN + ■ • • I (1) 

where the superscript refers to the number of derivatives or pion masses (chiral 
dimension) and the ellipsis denotes terms of chiral order four or higher. 
At lowest order, the Lagrangian in its relativistic form reads 

^i'l = * (*7^^M - ^'^A^ + Y^''^5"p) * (2) 
with 

D^ = d^, + F^ (3) 

r^ = ^(e^a^e + ^d,e) « ^ T . (tt X 9^7r) + . . . (4) 

u^. = ^{e^^.^ - ea^e^) « -^ T • 9^7r + . . . (5) 

/ti- 

U = ^i (6) 

For the parameters that occur in the first order Lagrangian, we use M^ = 
938.9187 MeV, U = 92.4 MeV, and gA = g-nNN Ul^'i-N = 1-29; the latter is 
equivalent to g^jq-^lA-ix = 13.67. 

Wp-will apply the heavy baryon (HB) formulation of chiral perturbation 
theoryliJ in which the relativistic Lagrangian is subjected to an expansion in 
terms of powers of I/Mat (kind of a nonrelativistic expansion), the lowest 
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In the relativistic formulation, the field operators representing nucleons, \E', 
contain four-component Dirac spinors; while in the HB version, the field oper- 
ators, N, contain Pauli spinors; in addition, all nucleon field operators contain 
Pauli spinors describing the isospin of the nucleon. 
At second order, the relativistic Lagrangian reads 

4 
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ttN 
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The various operators O^^' are given in Ref.E3. The fundamental rule by which 
this Lagrangian — as well as all the other ones — are assembled is that they must 
contain all terms consistent with chiral symmetry and Lorentz invariance 
(apart from the other trivial symmetries) at a given chiral dimension (here: 
order two). The parameters c; are known as low-enery constants (LECs) and 
must be determined empirically fron^fits to vrA^ data. We use the values 
determined by Biittiker and Mei6ncr,li3 which are (in units of GeV~^), 



ci = -0.81, 



C3 = -4.70 , 



C4 = 3.40 ; 



(11) 



C2 will not be needed. 

The HB projected Lagrangian at order two is most conveniently broken 
up into two pieces. 
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Note that C 
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TrAf.fix 



is created entirely from the HB expansion of the relativistic 

(2) 



£^^ and thus has no free parameters ( "fixed" ) , while £^^ ^^ is dominated by 
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the new ttN contact terms proportional to the q parameters, besides some 
small 1/Mjv corrections. 

At third order, the relativistic Lagrangian can be formally written as 

23 

4'ir = E^'*of)vi/, (15) 

with the operators, 0\ , listed in Refs.tj'E3; not all 23 terms are relevant to 
our problem. Similar to the order two case, the HB projected Lagrangian at 
order three is, 

f(3) _ f(3) , A3) (-.(.-^ 

'-ttN ~ '-ttN, fix ^ -^ttAT, ct' y-") 



with 



£i%,fix^+^Nh-i<J-u)DN 



9A ^ 



8Mf, 



(a- D){u- D)+h.c. 



N 



9A ^ 



l&Ml 



[a -u) D"^ + h.c. N + ... , (17) 



the ellipsis standing for relativistic correction terms not needed here, and 
^ttN ct given in Refs.Oll3; the latter reference contains also a determination 
of the di LECs. 

3 Pion-Exchange Diagrams and Power Counting 

The nN Lagrangian constructed in the previous section is the crucial ingre- 
dient for the evaluation of the pion-exchange contributions to the NN in- 
teraction. Since we are dealing here with a low-energy effective theory, it is 
appropriate to analyze the contributions in terms of powers of small momenta: 
{Q/A^^Y, where Q is a generic momentum or a pion mass and A^^ « 1 GeV 
is the chiral symmetry breaking scale. This procedure has become known as 
Power Counting. For the pion-exchange diagrams relevant to our problem, 
the power u oi & diagram is determined by the simple formula 

v = 2x loops + ^(dj - 1) , (18) 

j 

where ^loops'' denotes the number of loops in the diagram, dj the number of 
derivatives involved in vertex j , and the sum runs over all vertices. 

The most important irreducible one-pion exchange (OPE) and two-pion 
exchange (TPE) contributions to the NN interaction up to order Q^ are shown 
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Figure 1. The most important irreducible one- and two-pion exchange contributions to the 
NN interaction up to order Q'^. Vertices denoted by small dots are from C L, while large 
dots refer to C L .. 

Triv , ct 



in Fig. 1; they have been evaluated by Kaiser et alB using covariant pertur- 
bation theory and dimensional regularization. In addition to the diagrams 
displayed, we take into account the relativistic corrections up to order three 
as implied by £^^ g^, Eq. (jl^) and £^j(, j^^, Eq. (p^), where the latter con- 
tributes only to the one-pion exchange — to the order we are working at. 

One- and two-pion exchanges are known to describe NN scattering in 
peripheral partial waves. Therefore, we show in Fig. 2 (solid line) predictions 
by the chiral model displayed in Fig. 1 (plus realtivistic corrections up to 
order three) for the phase shifts in G waves. To provide a comparison with 
conventional mesorj_theory, we show also the predictions for tt -I- 27r exchange 
by the Bonn modeltO (dashed line). The dotted line represents pure one-pion 
exchange. Note that our calculations with 27r models always include also the 
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Figure 2. G-wave phase shifts. The predictions from the chiral n>»|iel displayed in Fig. 1 
are shown by the solid curve and the ones from the Bonn modeltJ by the dashed curTpi 
The dotted curve is OPE. Solid dots rep^psent the Nijmegen multi-energy np analysisE3 
and open circles the VPI/GWU analysisjij 



iterated one-pion exchange. From Fig. 2 we can conclude that, in G waves 
(orbital angular momentum L = 4), there is good agreement between the 
chiral and conventional I-k model as well as the empirical phase shifts. This 
is also true for all partial waves with i > 4. 

The agreement deteriorates when proceeding to lower L. While in F 
waves the agreement between the chiral model and the empirical phase shifts 
is still fair, substantial discrepancies emerge in D waves as demonstrated in 
Fig. 3: the chiral 27r exchange is far too attractive — a fact that has been 
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Figure 3. Z)-wave phase shifts. Notation as in Fig. 2. 



noticed before.&Ll 

4 Short-Range/Contact Contributions 

To control the D (and lower) partial waves, we need (repulsive) short-range 
contributions. In the conventional meson model, these are created by the ex- 
change of heavy mesons (notably, the to meson). In chiral perturbation theory 
(xPT), heavy mesons have no place and the short-range force is parametrized 
in terms of contact potentials, which are organized by powers of Q. If Q 
is, e. g., a momentum transfer, i. e., Q = p ' — p, where p and p ' are the 
CM nucleon momenta before and after scattering, respectively, and 9 is the 
scattering angle, then, for even v, 



Q" r^icosdy^ with 



m < — . 
- 2 



Partial-wave decomposition for orbital-angular momentum L yields, 

Q''PLicose)dcose^O for L<-, 
1 2 



(19) 



(20) 



where Pj, is a Legendre polynominal. The conclusion is that for non-vanishing 
contributions in D waves {L — 2), i^ = 4 is required. Based upon invariance 
considerations, there are a total of 24 contact terms up to order Q^, which we 
all include in our model. The parameters of these terms have to be natural, but 
are otherwise not restricted and, thus, represent essentially free parameters. 
The ideas of xPT may suggest that, if contacts are included up to order 
four, then also the 27r contribution should be calculated up to order four. (We 
went up to third order in the previous section, cf. Fig. 1.) We have looked 
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into this issue and found that the current status of the chiral 27r exchange 
at order four is a Pandora's Box. Some contributions have been calculated, 
like the football diagram with both vertices taken from C\}^ ^.j, i. e., both 
vertices proportional to one of the LECs Ci . It turns outE3 that the attraction 
generated by this diagram alone is so huge that it essentially doubles the size 
of the chiral 27r exchange in peripheral partial waves, leading to very serious 
(irreparable?) discrepancies with the empirical phase shifts and the predic- 
tions by conventional meson-exchange models. Other contributions at order 
four are obtained by inserting for the solid-dot vertices in Fig. 1 interactions 
from £^j(, ^^ (proportional to the LECs di). Consistency requires that these 
diagrams are calculated together with corresponding two-loop contributions 
of order four. All this is very involved and it will create horrific mathemat- 
ical expression to just describe the long-range NN interaction. At this time, 
nothing is known about these contributions, but the most optimistic prognosis 
would be that the latter contributions will compensate the enormous attrac- 
tion from the cf football. But even if this optimistic scenario were to come 
true, the chiral model would no longer be practical for nuclear physics pur- 
poses. It may then be much more reasonable to use the 'resonance saturation' 
argument and return to the conventional meson models oLthe past, particu- 
larly, the beautiful and simple one-boson-exchange modehcj in this model the 
one-sigma exchange (the mathematical expression for which can be written in 
just one line) describes the entire 2tt exchange — quantitatively! 

In conclusion: for the time being, the only realistic avenue towards a 
quantitative NN potential based upon xPT is a 'split approach': chiral 27r up 
to order three, contacts up to order four. This is our model. 

5 The NN Potential 

Since the 27r exchange diagrams. Fig. 1, are calculated using covariant per- 
turbatioii-.theory,El it is appropriate to start from the Bethe-Salpeter (BS) 
equationEa which reads in operator notation 

T = V + VgT (21) 

with T the invariant amplitude for the two-nuclcon scattering process, V the 
sum of all connected two-particle irreducible diagrams, and G the relativis- 
tic two-nucleon propagator. The BS equation is equivalent to a set of two 
equations: 

T = V + VgT (22) 

V = V + V{g-g)V (23) 
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« V + VopE {G - g)VoPE (24) 

where the last hne states the approximation we are using, exhibiting the way 
we treat the 27r box diagram (Vqpe is the relativistic one-pion exchange and 
V contains aU the irreducible diagrams of the type displayed in Fig. 1, includ- 
ing relativistic corrections up to order three). This treatment avoids double 
counting when V is iterated in the scattering equation and is also consistent 
with the calculations of Ref.El. For the relativistic three-dimensionaJ propa- 
gator g, we choose the one proposed by Blankenbecler and Sugara (BbS)° 
which has the great practical advantage that the OPE (and the entire poten- 
tial) becomes energy-independent. Thus, we do not need the rather elaborate 
formalism of unitary transformations^ to generate energy-independence of the 
potential. 

Our full chiral NN potential V is defined by 

j-r/^i . f sum of irreducible 1 ^ ^ ,._. 

^b ' P) ^ I ^ + 27r contributions | + ^°"*^^*^ ' ^^5) 

where the first term on the r.h.s. is given by Eq. (p4). 

This potential satisfies the relativistic BbS equation, Eq. (p2[). If we define 
now. 



E. 



V 




V{p',p) ^ ,r-^V{p\p) W^ « ( 1 - ^-rrrf) V{p\p) (26) 



with Ep = \/ M'^ -V jP' , then V satisfies the usual, nonrelativistic Lippmann- 
Schwinger (LS) equation, 

t M 

T{p\p) = V{p',p) + / dV V{p',p ") jj^— T{p ",p) . (27) 

J p'' — p" + It 



Note that the correction term — (p'^ -(-p^)/4M^ in Eq. pfl) is included only 
for OPE; for TPE it would create contributions beyond the order to which we 
calculate and for contacts it creates either existing terms or goes beyond our 
accuracy. 

In summary, our chiral NN potential V is defined by Eq. (pq) with V as 
given in Eq. (Eq). Since V satisfies Eq. (|27|), it is suitable for application in 
conventional, nonrelativistic nuclear structure physics. 

Iteration of V in the LS equation requires cutting V off for high momenta 
to avoid infinities. Therefore, we regularize V in the following way: 

V{jS\p) ^^ V{p',p) e-(p'/^)" e-(P/^)'" (28) 

"^For a derivation of the BbS approach, see appendix A.l of Bef.MI. 
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V{p',p){l- 



2n 



(29) 



where the last equation gives an indication of the fact that the exponential 
cutoff does not affect the order to which we are calculating, but introduces 
contributions beyond that order. For the contact potentials, we use cutoff 
masses A which are partial wave dependent. One can show that in doing so 
we just generalize the above regularization concept in the following sense: 



V{p',p)^V{p',p){l~h 



\A 



2n 



(30) 



with 6, of 0(1). 

For OPE, we use n = 4 and A = 0.6 GeV; for TPE, n = 2 and A = 0.46 
GeV; and for the contact potentials, n — 2 (except for contacts of order Q^ 
where n — i) and A « 0.4 — 0.5 GeV. The number of cutoff parameters 
is 22 which, together with the 24 contact parameters, results in a total of 
46 parameters for our model. At first glance, thiH may sound a lot. Note, 
however, that the Nijmegen phasp shift analysial3 uses 40 parameters and 
that the high-precision potentialsEJ^EJEHl developed in the 1990's have between 
40 and 50 parameters. Thus, a precise fit of the NN data requires around 50 
parameters, unless a model has more predictive power than the meson model. 
An important comment that has to be made about the chiral NN potential 
is that it has very little predictive power — less than the meson model. In the 
light of this fact, the number of 46 parameters is no surprise. 



6 Results for the Two-Nucleon System 

6.1 Two-Nucleon Scattering 

In Figs. 4 and 5, we show the phase shifts of neutron-proton {np) scattering 
for lab. energies below 300 MeV and partial waves with J < 2. The solid line 
represents the result from the chiral NN potential developed in the present 
work. The reproduction of the empirical phase shifts by our model is excellent. 
For comparison, we also show the phasa shift predictions by two chiral models 
recently developed by Epelbaum et aln (dotted and dashed curves in Figs. 4 
and 5). In the upper part of Table 1, we give our results for the effective range 
parameters of the S waves which agree accurately with the empirical values. 
We note that our present chiral potential is charge-independent and adjusted 
to the np data. 
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Figure 4. Phase shifts for J < 1. The solid Une is the result from our chiral NN potential, 
while the dotted and dashed lines are the prediotions by two chiral models developed by 
Epelbaum et al. (NLO and NNLO, respectively) .y The notation for the empirical points is 
the same as in Fig. 2. 
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Figure 5. Mixing angle ei and J = 2 piiase parameters. Notation as in Fig. 4. 



6.2 The Deuteron 

The reproduction of the deuteron parameters is shown in the middle part 
of Table 1. We present results for two versions of our chiral NN potential, 
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Table 1. Two- and three- nucleon low-energy data. 





Idaho-A" 


Idaho-B" 


CD-Bonn 


n 


AVl|] 


Empirical'' 


Low-energy np scattering 










""^50 scattering length (fm) 


-23.75 


-23.75 


-23.74 -23.73 


-23.74(2) 


^Sq effective range (fm) 


2.70 


2.70 


2.67 2.70 


2.77(5) 


^Si scattering length (fm) 


5.417 


5.417 


5.420 5.419 


5.419(7) 


^5i effective range (fm) 


1.750 


1.750 


1.751 1.753 


1.753(8) 


Deuteron properties 










Binding energy (MeV) 


2.224575 


2.224575 


2.224575 2.224575 2.224575(9) 


Asympt. S state (fm-l/^) 


0.8846 


0.8846 


0.8846 0.8850 


0.8846(9) 


Asympt. D/ S state 


0.0256 


0.0255 


0.0256 0.0250 


0.0256(4) 


Deuteron radius (fm) 


1.9756^= 


1.9758= 


1.970= 1.971= 


1.9754(9)'' 


Quadrupole moment (fm^) 


0.281"= 


0.284= 


0.280= 0.280= 


0.2859(3) 


_D-state probability (%) 


4.17 


4.94 


4.85 5.76 




Triton binding (MeV) 


8.14 


8.02 


8.00 7.62 


8.48 



"Chiral NN potential of the present work. 

''For references concerning the empirical data, see Tables XIV andi^VIII of Ref. 
=With mesp«-exchange current (MEC) and relativistic corrections .r1 
''ReferenceEj. i— i 

=Including MEC and relativistic corrections in the amount of 0.010 fm^.Ej 



.0. 



dubbed 'Idaho-A' and 'Idaho-B'.HThe main difference between the two models 
is in the D-state probabiUty of the deuteron, Pp. Even though Po is not 
an observable, it is of theoretical interest since the binding energies of few- 
and many-nucleon systems depend on it (cf. triton results at the bottom of 
Table 1). As mentioned before, the predictive power of the chiral model is 
very limited and it is possible to construct chiral potentials that fit the ^Si, 
^Di, and ei phase parameters up to 300 MeV and the empirical deuteron 
properties accurately, but have D-state probabilities that range from 3 to 6%. 
Such a large variation of Pd is not possible within the nrcson model of nuclear 
forces. 

Remarkable are the results produced by our chiral potentials for the 
deuteron radius which agree accurately j-sfith the latest empirical value ob- 
tained by using the isotope-shift methodBa All NN potentials p£the past (Ta- 
ble 1 includes two representative examples, namely, CD-Bonnt3 and AiOStJ) 
fail to reproduce this very precise new value for the deuteron radius.EJ Our 
chiral NN potentials are the first to predict this value right. 

In Fig. 6, we display the deuteron wave functions derived from our chiral 



''We note that the phase shifts represented by the solid line in Figs 4 and 5 are for Idaho-B; 
however, the ones for Idaho-A are so close to Idaho-B that they could not be distinguished 
on the scale of the figure. 
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r(fm) 

Figure 6. Deuteron wave functions: the larger curves are 5-waves, the smaller ones D- 
waves. The solid line represents the wave functions derived from our chiral NN potential 
(Idaho-BVr-iThe daahad, dash-dojrt«d, and dotted lines refer to the wave functions of the 
CD-Bonnll3, Nijm-llia, and AVlstil potentials, respectively. 



potential (Idaho-B) by the solid line and compare them with wave functions 
based upon conventional NN potentials from the recent past. Characteristic 
differences are noticeable; in particular, the chiral wave functions are shifted 
towards larger r which explains the larger deuteron radius. 

Concerning the triton binding energy predictions given at the bottom of 
Table 1, we like to comment that the results for Idaho- A and B are obtained in 
a 34-channel Faddeev calculation with no charge-dependence (i. e., using the 
np potential throughout), while the corresponding calculations with CD-Bonn 
and AV18 take charge-denpendence into account. 

7 Summary and Conclusions 

We have constructed an accurate chiral NN potential. The model includes 
one- and two-pion exchange contributions up to chiral order three and contact 
terms (which represent the short range force) up to order four. Within this 
framework, the NN phase shifts below 300 MeV lab. energy and the properties 
of the deuteron are reproduced with high-precision. 
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Due to the very quantitative nature of this new chiral NN potential, it rep- 
resents a rehable and promising starting point for exact few-body calculations 
and microscopic nuclear many-body theory. 
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